Evolution of reference networks with aging 
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We study the growth of a reference network with aging of sites defined in the following way. 
Each new site of the network is connected to some old site with probability proportional (i) to the 
connectivity of the old site as in the Barabasi- Albert's model and (ii) to t~", where r is the age of 
the old site. We consider a of any sign although reasonable values are < a < oo. We find both 
from simulation and analytically that the network shows scaling behavior only in the region a < 1. 
When a increases from —oo to 0, the exponent 7 of the distribution of connectivities {P{k) oc k^'' for 
large k) grows from 2 to the value for the network without aging, i.e. to 3 for the Barabasi- Albert's 
model. The following increase of a to 1 makes 7 to grow to cxd. For a > 1 the distribution P{k) is 
exponentional, and the network has a chain structure. 
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The explosion of the general interest on the problem 
of the structure and evolution of most different networks 
is connected not only with the sudden understand- 
ing that our world is in fact a huge set of various networks 
(with the most striking exarnples of Web and neural net- 
works - see the papers [|[|[||,^ and references therein) 
but also with the recent finding that many networks obey 
scaling laws_[0,|j. That was the reason to renew old 
studies p0|-|13[. One may note also that the growth of 
networks is only a particular kind of fascinating growth 
processes [p^-p^. 

The simplest kind of growing networks is the reference 
networks |1^ in which new links appear only between 
new sites and the old ones but they never appear be- 
tween old sites. Therefore, in the reference networks, the 
average connectivity of old sites is always higher than 
that one of younger sites. The well-known example of 
the reference networks is the network of citations of sci- 
entific papers in which each paper is a site of the cor- 
responding net and links are the references to the cited 
papers. 

A clear beautiful model of a reference network 
that shows scaling behavior was recently proposed by 
Barabasi and Albert |^. In their network, each new site 
is connected with some old site with probability propor- 
tional to its connectivity k. In this case, the distribu- 
tion of the connectivities in the large network (i.e. one of 
the most considerable characteristics of the structure and 
evolution of networks) shows a power-law dependence 
P{k) oc k^"' with the exponent 7 = 3. 

However, in real reference networks aging usually oc- 
curs: alas, we rarely cite old papers! One may ask, how 
do the structure of the network change if the aging of 
sites is introduced, i.e. if the probability of connection of 
the new site with some old one is proportional not only 
to the connectivity of the old site but also to the power 
of its age, t~", for example? This question is quite rea- 
sonable: indeed. Fig. 1 demonstrates that the structure 
of the network depends obviously on a. We show be- 



low both numerically and analytically that this change 
is dramatic: the scaling disappears when a > 1, and the 
exponents of the scaling laws depend strongly on a in 
the range a < 1. This result could not be foreseen be- 
forehand: the scaling is very sensitive to changes of the 
model. For example, as it is noted in the scaling dis- 
appears if the probability of the connection to an old site 
is proportional to fc', e 7^ 1. 

Here, we consider exclusively the network, in which 
only one extra link appears when the new site is added, 
since the results for the exponents do not depend on the 
number of links which are added each time . 

Let us start from the simulation which turns to be easy 
for the problem under consideration because we study 
only the characteristics dependent on the connectivity. 
They are: (i) the distribution of connectivities in the 
network P[k,t) in the instant t (only one site is in the 
network at i = 0, and one more site is added in each in- 
stant) and (ii) the mean connectivity k(s, t) of the site s 
(0 < s < i) in the instant t, i. e. the local density of the 
connectivities. If one is interested only in these quanti- 
ties, there is no need to keep a matrix of connections in 
memory, and the simulation is very fast. 

The results of the numerics are shown in Figs. 2-5. 
Although only the region a > seems to be of real 
significance, we consider also negative values of a since 
they do not lead to any contradiction. One may see 
clearly from the figures that P{k, t) oc k^^ for large k 
and k{s, t) oc for small s, where (3 is the other scal- 
ing exponent, only for a < 1. As it should be, we get 
the values 7 = 3 and /3 = 1/2 for a = Note that, 
at < a < 1, the deviations of the dependence logP(fc) 
vs. log k for small k are stronger as in the real reference 
network ||^, than those ones at a = 0. At a > 1, P{k) 
turns to be exponential, and the mean connectivity tends 
to be constant at large s. Thus, while a changes from 
to 1, 7 grows from 3 to 00, and (3 decreases from 1/2 to 
0. One sees from Fig. 3, that, for positive a the log P{k) 
vs. log k dependence looks more curved than those ones 
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for a = like in the real reference network 

The simulations demonstrate also that 7 has a ten- 
dency to decrease from 3 to 2, and /3 to grow from 1/2 
to 1 when a decreases from to —00. 

Let us show now how these results may be described. 
We use an effective medium approach, inspired by 
We shall explain why it gives so good results for the ex- 
ponents later but first we outline the applied theory for 
the problem under consideration. 

The main ansatz of the effective medium approach in 
our case is the approximation of the probability P{k, s, t) 
that the connectivity of the site s at the moment t is equal 
to k by the ^-function: 



P{k,s,t) = 6{k - k{s,t)). 



(1) 



Thus, the exact function P{k,s,t) is assumed to be 
enough sharp. We suppose also that we may use a con- 
tinuous approximation. These two strong assumptions 
lead immediately to the following equation: 



dk{s,t) 
dt 



k{s,t)(t-s)-^ 
t _ 

J duk{u,t){t - u)-" 




k{t,t) = 1. 



(2) 



Here, the boundary condition means that only one link 
is added each time. 

One may check that Eq. is consistent. Let us apply 

ds to it. Then, 



^^J-kis,t) 



d 



dt dt JO 



dsk{s,t) - k{t,t) ^ 1, (3) 



and we get immediately the proper relation 



dsk{s,t) = 2t, 



(4) 



i. e. the sum of connectivities equals the doubled number 
of links in the network. 

We search for the solution of Eq. ^ in the scaling 
form 



k{s,t) = K{s/t) , s/t = 

which is consistent also with Eq. (H). Then Eq. 
becomes 



(5) 



1 -,-1 

dC/^(C)(i-C)-" 



(6) 



where /3 is a constant, which is unknown yet. We 
shall see soon that this constant is the exponent of the 
mean connectivity. Eqs. and give the relation 
J^dCK{C)=2. 

The solution of Eq. (||) is 



k(^) — B exp 



d^ 



C(i-e)" 



(7) 



where _B is a constant. The indefinite integral in Eq. (j^) 
may be taken: 
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fc=0 



fc!(fc + 1): 



;a{a + l)...{a + k)£, 



k+l _ 



lne + a3-F2(l,l,l + a;2,2;0, 



(8) 



where 3-F2( ,,;,;) is the hypergeometric function [ p^ . 
Recalling the boundary condition k{1) = 1, we find the 
constant B. Thus the solution is 

«(0 = 

^_^(C+V,(l-o))^-/3 ^^^(1^ 1^1 + ^. 2^ 2; ^)] , (9) 

where C = 0.5772 ... is the Euler's constant and '4){ ) is 
the T/i-function. Now we see that the constant (3 indeed 
is the exponent of mean connectivity, since k(^) ^ f if 
^ — > 0. The transcendental equation for /3 may be writ- 
ten if one substitutes Eq. (0) into the right side of Eq. 



^-1 ^ g-/3(C-(-^(l-a)) X 



C^(i-C)" 



exp [-/3aC 3^^2(1, 1,1 + a; 2, 2; 0], (10) 



that is our main equation. 

Before we shall find the solution of Eq. (|l0|), let us 
show how the exponents (3 and 7 are related if one applies 
Eq. (|l|) of the effective medium approach. k{s,t) oc s^^, 
so s oc k^'^/f^ . Hence, k~'^ oc P{k,t) oc ds/dk oc k^^^'^/f^^ 
and one obtains that 



7=1 + l/P. 



(11) 



The solution of Eq. ( p^ ) exists in the range — 00 < 
a < 1. The results of the numerical solution are shown 
in Fig. 5. One may also find /3(ck) and 7(a) at a — > 0: 

/3 = i - (1 -ln2)a , 7 ^ 3 + 4(1 -ln2)a, (12) 

where the numerical values of the coefficients are 1 — 
ln2 = 0.3069... and 4(1 -In 2) = 1.2274.... We used the 

relation 3i^2(l, 1'1;2J;C) = i*2(C)/C = CVfc')/C 
while deriving Eq. (|l2|) . Here, Li2{ ) is the polylogarithm 
function of order 2 [ |18[ . 

In the limit of a 1, using the relation 

3^2(1, 1, 2; 2, 2; C) = - ln(l - C)/C, we find 



(} = ci(l - a) , 7 = 



1 



ci 1 — a 



(13) 
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Here, a = 0.8065 . . . , = 1.2400 . . .: the constant ci 
is the sohition of the equation 1 + 1/ci = exp(ci). Note 
also that (3 1 and 7 — > 2 in the hmit a —>■ —00. 

One sees from Fig. 5 that the resuhs of the simu- 
lation and of the analytical calculations are in qualita- 
tive correspondence. Deviations may be noticed only 
in regions where simulations can not provide sufficient 
precision. Why does such a coarse theory demonstrate 
so close agreement with the simulation? One may show 
p9[ that the function P{k,s,t) of k is functionally very 
sharp as compared with the studied long-tailed distribu- 
tion P{k,s). Therefore, the main anzats Eq. (|l|) gives 
excellent results. In fact, Eqs. (0) and (||) provides us 
easily with all known results on reference networks. 

One may imagine from Fig. 1 that the point a ^ 1 
(/3 ^ and 7 00) at which the scaling collapses 
marks the transition from the multidimensional network 
for a < 1 to the chain structure. Fig. 1 demonstrates also 
that, in the case of a ^ —00 {f3 ^ 1 and 7 — > 2) all sites 
of the network are connected with the oldest one. The 
behavior of the considered quantities near these points 
evokes associations with the lower and higher critical di- 
mensions in the theory of usual phase transitions pO|] . 

There are some possibilities to change the exponents 
of the network without aging One may check that, 

in these cases, the range of variation of the exponents /3 
and 7 when a changes from — cxd to 1 is the same as in 
the present Letter. 

In summary, we have considered the reference network 
with the power-law (t~") aging of sites. We have found 
both from our simulations and using the effective medium 
approach that the network shows scaling behavior only 
in the region a < 1. We have calculated the expo- 
nents of the power-law dependences of the distribution 
of connectivities P{k,t) oc k^'^ and of mean connectivity 
k{s,t) oc and have shown that they depend crucially 
on a. 

The following questions remain open. Is our unproved 
idea about the nature of the threshold point a = 1 rea- 
sonable? Do the analogies with the lower and higher 
critical dimensions exist indeed? What other quantities 
of the network demonstrates the scaling behavior? 
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a = -10.0 a= 0.0 




a = 0.5 a = 1 .0 




a = 2.0 a = 10.0 

FIG. 1. Change of the network structure with increase 
of the aging exponent a. The aging is proportional to t~" , 
where r is the age of the site. The network grows clockwise 
starting from the site below on the left. Each time one new 
site with one link is added. 



FIG. 2. Mean connectivity vs. number s of the site for 
several values of the aging exponent a: 1) a = 0, 2) a = 0.25, 
3) a = 0.5, 4) a = 2.0. The network size is t = 10 000. 
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FIG. 3. Distribution of the connectivities for several val- 
ues of the aging exponent: 1) a = 0, 2) a = 0.25, 3) a — 0.5, 
4) a — 0.75. The inset shows logP(fc) vs. k for a = 1,5, and 
10. Note that, in the later case, all three curves are nearly 
the same. 
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FIG. 5. Exponent 7 of the connectivity distribution vs. 
aging exponent a. The points show the results of the sim- 
ulations. The line is the solution of Eq. (^o|) with account 
for Eq. (|ll|). The inset depicts the analytical solution in the 
range —5 < a < 1. 



